We present the basic features and preliminary results of the interface between our spectro-photometric model GRASIL (that calculates galactic SED from the UV to the sub-mm with a detailed computation of dust extinction and thermal reemission) with the semi-analytical galaxy formation model GALFORM (that computes galaxy formation and evolution in the hierarchical scenario, providing the star formation history as an input to our model). With these two models we are able to synthetize simulated samples of a few thousands galaxies suited for statistical studies of galaxy properties to investigate on galaxy formation and evolution. We find good agreement with the available data of SED and luminosity functions.